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Abstract 

A large class C of relativistic diffusions with values in the phase-space of special relativity was 
introduced in [5] in order to answer some open questions concerning the asymptotic behaviour 
of two examples of such processes [IT] [20] [21] . In particular, the equilibrium measures of these 
diffusions were explicitly computed, and their hydrodynamic limit was shown to be Brownian. 
In this paper, we address the question of the trends to equilibrium of the diffusions of the whole 
class C. We show the existence of a spectral gap using the method introduced in [7] and deduce 
the exponential decay of the distance to equilibrium in L 2 — norm and in total variation. A 
similar result was obtained recently in [12] for a particular process of the class C. 

Keywords: Rclativisic diffusions, Relativistic Ornstein-Uhlcnbcck process, Equilibrium measure, 
Spectral gap, Lyapounov function. 
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1 Introduction 

The study of stochastic processes in the framework of special relativity goes back to the 1960s and 
the pionneering work of Dudley [T8][T!|]. Since the late 1990s, there has been a renewed interest in 
the subject with the work of Debbasch and his co-authors on the Relativistic Ornstein-Uhlenbeck 
Process [9] [TOJ [T6] [17] , then those of Dunkel and Hanggi on the so-called "relativistic Brownian 
motion" [30] [H]. The notion of relativistic diffusion has been extented to the realm of general 
relativity in [33] [T5] and the literature on the topic is now thriving, both in Mathematics, see for 
example [33] H [3J and in Physics [Ml E] etc. 



In this article, we address the question of the trends to equilibrium for a large class C of relativis- 
tic diffusions with values in the phase-space of special relativity, i.e. the unitary tangent bundle of 
Minkowski space-time. This class of processes was introduced in [5] in order to answer several open 
questions concerning the long-time asymptotic behavior of two examples of such diffusions, namely 
the ones considered in [!()] [201 EI] ■ To our knowledge, the class C includes most of the Minkowskian 
diffusions introduced in the physical literature. 
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We show that under the same mild hypotheses as in [5], and for all diffusions of the class C, the 
equilibrium measure of the "momentum subdiffusion" satisfies a Poincare inequality and we thus 
deduce that the rate of convergence to equilibrium is exponential both in L 2 — norm and in total 
variation. A similar result was obtained recently in [12] for a particular process of the class C: the 
relativistic diffusion process associated to the kinetic relativistic Fokker-Planck equation considered 
in [TJ [SO] [HJ [51] . The method we follow here is the one developped in [7] [5] , which generalize the 
classical Bakry-Emery criterion when the potential associated to the equilibrium measure is not 
strictly convex. It is based on the existence of a Lyapounov function associated to the infinitesimal 
generator of the diffusion. 

The structure of the article is the following: in the next section, we introduce some notations 
and we recall the definition of the class C of relativistic diffusions considered in the sequel. In 
Section [3] we state our results concerning the trends to equilibrium of the momentum components 
of the diffusions. The last section [4] is devoted to the proof of our main result, namely the Poincare 
inequality satisfied by the equilibrium measure. 

2 The class C of relativistic diffusions 

Fix d > 1 an integer, and denote by ||x|| = \J\x x | 2 + . . . + |.T d | 2 the Euclidian norm of a vector 
x € K d . Let M. 1,d denote the Minkowski space of special relativity. In its canonical basis, denote 
by x = (x^) = (x°, x l ) = (x°, x) the coordinates of the generic point, with greek indices running 
0, .., d and latin indices running 1, .., d. The Minkowskian pseudo- metric is thus given by 

d 

ds 2 = \dx°\ 2 -J2\dx l \ 2 - 

i=l 

The world line of a particle with positive mass m is a timclikc path in which we can always 
parametrize by its arc- length, or proper time s. So the moves of such particle arc described by a 
path s i— > (x^) in K 1,d , having momentum p = (p s ) given by p = (p^) = (p°,p l ) = (p°,p), where 
d-Ps '■= fndx^/ds, and satisfying the pseudo-norm relation 

|p°| 2 -||p|| 2 = m 2 . 

We shall consider here future directed world lines of type (t, x t ) t >o, and take m = 1. Introducing the 
velocity v = (v 1 , . . . , v d ) by setting v l := dx % /dt, and working with the usual spherical coordinates 
(r,6) el + x S^" 1 , r := ||p|| and 6 := p/r =: (0 1 , . . . , d ) , we get at once: 

| p° =£ = v / IT^=(i-|v|T 1/2 , 

{ p = \Jl + r 2 (l,v). 

Thus, a full space-time trajectory (xt,pt)t>a = (t, ^-t,Pt, Pt)t>o which takes values in the positive 
part of the unitary tangent bundle T+M}' d , is determined by the mere knowledge of its spacial 
components (xt,pt)- We can therefore, from now on, focus on spacial trajectories t i— > (xj,pt) 
which take values in the Euclidian product R d x R d . 
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2.1 Definition of the class C 

Let us first recall the definition of the class C introduced in [5j . 

Definition 1. The relativistic diffusions of the class C are the processes of type (x t ,pt)t>o = 
(i,x t ,j»(, pt)t>o i n TIM. 1 ' , where the associated spatial process (x t ,p 4 ) t >o is itself a diffusion, so- 
lution of a stochastic differential system of the form, for 1 < i < d: 

( dx\ = f(n)p\dt 

* { dpi = -6(rOpjdt + CT(rO(/3[l+7 7 (r t ) 2 ])~ 1/2 [d^ t < +j ? (r t )ejd«; t ] 

where the real functions f , b, a , r] are continuous on and satisfy the following hypotheses, for 
some fixed e > 0: 



2rb(r) 

(n) { " ' »v -aHr) 

lim e 6 r f(r) = for some e' < 



a > e on K + ; g(r) := — > e for large r ; 



r— >oo 



In the definition above, W := (W 1 , . . . , W d ) denotes a standard d-dimensional Euclidian Brow- 
nian motion, w denotes a standard real Brownian motion, independent of W, and /3 > is an 
inverse heat parameter. 

Example 1. In the simplest case of constant functions f , b, a, and rj = 0, the process (x t ) t >o is an 
integrated Ornstein- Uhlenbeck process. The process considered by Debbasch et al. in JPj 1 101 \16l 
they call Relativistic Ornstein- Uhlenbeck Process (ROUP), corresponds to: 

f(r) = b(r) = (1 + r 2 )- 1/2 , a(r) = y/2, n = 0, 

and the process considered by Dunkel and Hdnggi \2(A \ 21^ corresponds to: 

f{r) = (1 + r 2 )- 1 / 2 , b(r) = 1 - dp- 1 (1 + r 2 )" 1 / 2 , 



o~(r) = v 2Vl + r 2 , n(r) = r. 
2.2 Infinitesimal generator of the momentum diffusion 

If (xt,pt)t>o = (t, xt,p° , p*)t>o is a relativistic diffusions of the class C, then the process (pt)t>o is 
itself a diffusion process, we will call the momentum diffusion. In spherical coordinates p = (r, 9) € 
M+ x § d_1 , its infinitesimal generator is given by 

a 2 (r) 

C a 2 := C r + Asd-i, (1) 

where A§<j-i denotes the usual spherical Laplacian on S d_1 and C r is the infinitesimal generator of 
the radial process r t := ||pt||: 



2/3 



x i /\ 2 +/ g g( r ) 



r 1 + rj(r) 2 



0, 
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Let us use the same notations as in [5], that is: 

M(0 " ^ (/ p(l+t(p) 2 ) ' P ) ' G{r) ' = I 9{P)dPl 
and introduce the functions V : R+ -t K and U : R d — > R: 



V(r):= /^x f ^ dp + ^r), C/(p) := V(||p||). 
1 P 1 + 



Then, if A and V are the usual Laplacian and gradient in M. d , the generator C a 2 at p can be 
re-written under the familiar form: 

£<x2 := ^rjp^- x A with £ := A - Vt/(p).V . (2) 

In the sequel, we denote by T the operator "carre du champ" associated to £, that is for good 
functions / and g in the domain of C: 

T(f,g):=C(fg)-C(f)g-fC(g). (3) 

3 Statement of the results 

We can now state our results concerning the trends to equilibrium of the momentum diffusion. In 
the following lemma, we explicit the invariant (or equilibrium) measure of the process (pt)t>o- 

Lemma 1. Let (£, x t ,pj , pt)t>o & diffusion of the class C. Then, the momentum subdiffusion (pt)t>o 
with values in M. d is ergodic and its invariant probability measure is given by: 

1 e~ u ^ 

" ip) := z x ^» dp ' (4) 

where dp denote the Lebesgue measure in M. d and Z is a normalizing constant. 

Remark 1. If the Euclidian space R d is endowed with the usual spherical coordinates (r,6) £ 
R + x § d_1 , the equilibrium measure v reads 

1 e- v ^ , n 1 ., , e"^ G W 

K^, °) = 7 x —qj-t-t drdd = - x M ( r )«»-i __drdfl, (5) 
Z (T z (rj Z a z (r) 

where dO is the uniform measure on S d_1 . 

Example 2. In the case of the ROUP and the diffusion of Dunkel and Hanggi, the invariant 
measure is the Jiittner distribution which is the equivalent of the classical Maxwell measure in the 
framework of special relativity: 

„(p) = I x e-'Vi+IIPlPdp, i.e. v (r,0) = -er^^^r^drdO. (6) 

z z 
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Proof. Fix a smooth, bounded test function / : M. d — > R. One easily check that 

so that the measure v is invariant for the process (pt)t>o- Under the hypotheses {%), one has 
min(l,r) < fj,(r) < max(l,r), a(r) > e for all r > and G(r) > er/2 for r sufficiently large so that 

-£/(p) /• , „-^G(r 



J cr 2 (||p||) J cr 2 (r) 



The measure 1/ being finite, the process (pt)t>o is ergodic. □ 

The next theorem and corollary establish that the equilibrium measure v satisfies a Poincare 
inequality, so that the rate of convergence to equilibrium of the momentum process is exponential 
in L 2 (i/)— norm. Here and in the sequel, vf denote the integral of the function / against v. 

Theorem 1. There exists a positive constant C2 such that for all f € L 2 (V): 

varAf) ■•=\\f-yf\\h {v) <C2 J T(J,f)du, 

Remark 2. Following the same method as in the proof of Theorem]]] in Section^ below, it can 
be shown (see JTjj) that the Poincare inequality even holds in lJ-(y), namely there exists a positive 
constant c\ such that for all g 6 L 1 (y) with zero median, one has: 



\g\\u(v) < ci J \Vg\dv. 



From the above theorem, by the well known equivalence between Poincare inequality and expo- 
nential decay of the L 2 (^)— distance to equilibrium, see for example Theorem 2.5.5 of [2], we deduce 
the corollary: 

Corollary 1. Let (t, Xj,Pj, Pt)t>o be a diffusion of the class C and let Pt be the Markov semi-group 
associated to the momentum subdiffusion (pt)t>o- Then, there exists a positive constant c^ such 
that for all t > 0, and for all f £ L 2 (i/): 

\\Ptf-vf\\hw<e- t '°'\\f- V f\\l, (v y 

Moreover, as it is shown in [T3], from the exponential decay of the L 2 (j/)— distance to equilib- 
rium, one can derive the exponential decay of the distance to equilibrium in total variation. Namely, 
if P*v denotes the law of pt with initial distribution v, we have: 

Corollary 2. Suppose that the initial distribution of (pt)t>o can be written v = hdv, with h € 
L 2 (i/). Then, for all t > 0, the distance in total variation satifies 

\\P*V-v\\ TV = ||P;/i-l|| L i W < e-'/^p- 1|| L2M . 
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The proof of Theorem [T] is given in the next section. As told in the introduction, the method we 
follow is an adaptation of the method developped in [3 |8] . It consists in expliciting a Lyapounov 
function for the generator L a i and then use this function to derive a Poincare inequality for the 
equilibrium measure v. The slight change from the original proof comes from the fact that in their 
original paper [7], Bakry et al. consider a framework where there is only an additive noise, i.e. 
a = 1 and r\ = 0. For the sake of completeness, a self-contained proof of the theorem in our context 
is given below. 

Remark 3. The exponential decay in Corollary^ is not trivial, even in the simplest examples of 
diffusions belonging to the class C. For example, in the case of the ROUP, the process (pt)t>o * s 
solution of the stochastic differential equations system: 

d Pt = P * 

Vi + IIPtll 2 

When \ \pt \\ is large, the drift is only linear in t so that the amount of time needed to return to zero 
may be very large compared to the case of the usual Ornstein- Uhlenbeck process where the drift is 
linear in p t . Equivalently, from an analytic point of view, the potential 



[/( P ) = ni|p|l) = Vi + llp|| 2 

is not strictly convex, so that the classical methods such as the Bakry-Emery criterion do not apply. 
However, such a result is not surprising since the exponential decay in If 1 — norm is well known (see 
for the very similar one dimensional process X t solution of the stochastic differential equation 

dX t = -sign(X t )dt + dW t , 

where the associated potential x i— > \x\ is not stricly convex either. In the last decade, much progress 
has been made to weaken the assumptions under which a measure in Euclidean space satisfies a 
Poincare inequality. Among recent advances, the method developped in [?J W$ provides an elegant 
and efficient way to get strong results, for example it applies to general log— concave measures. 



4 Proof of the results 

4.1 Existence of a Lyapounov function for the generator £ a 2 

We now give the proof of the theorem [T] stated in the previous section. We first explicit a Lyapounov 
function associated to the infinitesimal generator £ a 2 given by the equation J2]). We denote by 
B(0, R) the Euclidian ball centered at the origin and of radius R in R d . 

Lemma 2. There exists a smooth function W : M. d — > M, and some constants a>0,7>0, R > 
such that for all p £ l d : 

1. W(p) > 1 ; 

2. |t5t(p)| is bounded ; 

3. 4^(p)<-^(p)+ 7 l BM (p). 
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Proof. Consider a smooth function W of the form e c " p " for ||p|| > R and such that IU(p) > 1 for 
all p s R d , where the two constants c > and R > will be fixed later. For ||p|| > R, one has 

C^W( P ) = « x c + c - W(p).J^) WC P). 

In our case, we have Vt/(p) = VV(||p||) = U'(||p||).p/||p|| so that 

Ce*W(p) = x C + c - V'(HpII)) W(p). 



2/3 

Under the hypotheses (H), <7 2 (||p||) > e 2 for all p, and for large enough ||p|| we have: 

miPii) = ^ x - f il jfl 2 + M\P\\) > 

HpII i + ^dlPlI) 

Thus, taking c sufficiently small and R large enough so that + c < 4p, we get for ||p|| > R: 

e 3 

£ a sW(p) < — aW(p), where a = c x — . 
Finally, for some non- negative constant 7, we have for all p £ R d : 

4^(P) < -«VF(p)+ 7 l B (o,i?)(p)- 

Moreover, by construction ^p-(p) is bounded. □ 

Remark 4. Since the Lyapounov function W satisfies W(p) > 1 for all p € R rf , i/ie point 3. in 
Lemma\^ can be re-written as: 

i caw) 7 1 1 - 

Therefore, since W{\\p\\) goes to infinity with ||p||, i/ie iwo ratios — C " 2 ^ V ^ and — are positive 
for ||p|| Zarge enough. 

4.2 Proof of the Poincare inequality 

We now give the proof of the Poincare inequality stated in Theorem [TJ In the sequel dX denote the 
Lebesgue measure in R d . 

Proof. If g is a smooth function in L 2 (i/), we have vaiv(g) < J(g — c) 2 dv for all real constants c. 
Let c be such a constant and define / := g — c. Using the inequality ([7]), we have 

a 1 /- £ ff 2(WQ j2j 7 f f 



A B 



The first term A could be infinite depending on the behavior of C a 2 (W)/W at infinity. Since we do 
not impose any integrability condition on this ratio, we have to restrict ourself in the calculation 
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below to the case where / = (g — c)x, where x is a smooth, non-negative, compactly supported 
function such that 1b(o.r) X < !• The general case is then obtained by taking a sequence of 
functions \n such that lB(o,nR) < Xn < 1, |V%„| < 1, and go to the limit, which is allowed thanks 
to the mononicity noticed at the end of Remark |4] Since the generator C is symmetric with respect 
to e~ u dX, we have 



A = 



CAW) 



W 



fdv = 



1 



2[3Z 



w 



-fe~ u dX 



1 



2/3Z 
1 

~2f3Z 
1 

2PZ 



( C- ) VWe~ U dX 
W 



yw 



•' v/ ir ' r ' /A 

\Vf\ 2 e~ u dX- [ V/-/ 



vw 



w 



w 



- u d\ 



~ u d\ 



and therefore 



A> — 



(9) 



Let us consider now the second term B in the right hand side of (HJ. It is well know that the 
measure v satisfies a Poincare inequality in the ball B(0,R), i.e. for a positive constant kr, we 
have: 

/ \ 2 

1 



/ dv < n R 



B(0,R) 



|V/| 2 ^+--— — I / fdv 

B(0,R) K)) \Jb{0,R) 



We choose c = v(B(0, R)) 1 x J B ^ Q R ^ gdv : so that the last term vanishes ; thus using the fact that 
\W\ > 1 and a 2 > e 2 , we have 



B < KR 



\Vf\ 2 dv< 



B(0,R) 



KR 

e 2 Z 



\Vf\ 2 e- u d\. 



(10) 



B(0,H) 



Putting ([9]) and (|T0l) together, we conclude that 



7 K fl 

2a/3Z '" ae 2 Z 



/ |V/| 2 e"^A 



1 



7 K i? 



2a/3Z ae 2 Z 

1 7 K fl 



\Vg\ 2 e~ u dX 




2afiZ ae 2 Z / 

1 7«fl 
2afiZ ae 2 Z / 



20™'* 
?{g,g)di> 



In other words, we have shown the desired Poincare inequality with the (very non optimal) Poincare 
constant 



C2 



Q 



□ 
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